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Abstract 

A new formulation of the thermodynamic field theory (TFT) is 
presented. In this new version, one of the basic restriction in the 
old theory, namely a closed-form solution for the thermodynamic field 
strength, has been removed. In addition, the general covariance prin- 
ciple is replaced by Prigogine's thermodynamic covariance principle 
(TCP). The introduction of TCP required the application of an appro- 
priate mathematical formalism, which has been referred to as the iso- 
entropic formalism. The validity of the Glansdorff-Prigogine Universal 
Criterion of Evolution, via geometrical arguments, is proven. A new 
set of thermodynamic field equations, able to determine the nonlinear 
corrections to the linear (" Onsager" ) transport coefficients, is also de- 
rived. The geometry of the thermodynamic space is non-Ricmannian 
tending to be Ricmannian for hight values of the entropy production. 
In this limit, we obtain again the same thermodynamic field equations 
found by the old theory. Applications of the theory, such as transport 
in magnetically confined plasmas, materials submitted to temperature 
and electric potential gradients or to unimolecular triangular chemical 
reactions can be found at references cited herein. 

1 Introduction 

It is known that, close to equilibrium, the transport equations (i.e. the flux- 
forces equations) of a thermodynamic system are provided by the Onsager 
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relations. Indicating with X^ and the thermodynamic forces and fluxes, 
respectively, the Onsager relations read 

Jfl = T 0flu X U (1) 

where t ^ u are the transport coefficients. We suppose that all quantities 
involved in Eqs (pQ) are written in dimensionless form. Matrix tq^ u can be 
decomposed into a sum of two matrices, one symmetric and the other skew- 
symmetric, which we denote with and foUis, respectively. The second 
principle of thermodynamics imposes that L^ u be a positive definite matrix. 
In this equation, as in the remainder of this paper, the Einstein summation 
convention on the repeated indexes is adopted. The most important prop- 
erty of Eqs (H|) is that near equilibrium, the coefficients t^ v are independent 
of the thermodynamic forces, so that 

|^ = (2) 

The region where Eqs ([2]) hold, is called Onsager's region or, the linear 
region. Many important theorems have been demonstrated for thermody- 
namic systems in the Onsager region. Among them, the most important 
one is the Minimum Entropy Production Theorem, showed by Prigogine in 
1945-1947 [l|. This theorem establishes that, in the Onsager region, for a — a 
processes^, a thermodynamic system relaxes towards a steady-state in such 
a way that the rate of the entropy production is negative 

da (da \ 

— < — = at the steady state (3) 

dt ~ \dt J J w 

where a = L^X^X" indicates the entropy production and t is time. In 
1954, Glansdorff and Prigogine demonstrated a moregeneral theorem, valid 
also when the system is out of the Onsager region [3j. They showed that, 
regardless of the type of processes, a thermodynamic system relaxes towards 
a steady-state in such a way that the following quantity V is negative 



V = J ^~ d ~f - [V = at the steady state) (4) 



1 Here, we adopt the De Groot-Mazur terminology Q|: when the velocity's distribution 
function of particles is an even (odd) function of the velocities of particles, the processes 
is referred to as a a-process (b-process) . It is possible to show that this definition implies 
that the a-processes only involve the symmetric part of the Onsager matrix whereas the 
6-processes only the skew-symmetric one. 
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Inequality Q reduces to inequality Q for a - o processes in the Onsager 
region. For spatially-extended systems, the expression in Eqs. (j3|) should be 
replaced by 



where dv is the volume element and the integration is over the entire space f2 
occupied by the system in question. J^{r,t) and X^{r,t) denote the space- 
time dependent fluxes and forces, respectively. The inequality expressed in 
@ [or in ([5|)] is referred to as the Universal Criterion of Evolution and it is 
the most general result obtained up to now in thermodynamics of irreversible 
processes. Out of the Onsager region, the transport coefficients may depend 
on the thermodynamic forces and Eqs ([2]) may loose their validity. Trans- 
port, in the nonlinear region, has been largely studied, both experimentally 
and theoretically. In particular, many theories, based on the Fourier ex- 
pansion of the transport coefficients in terms of the thermodynamic forces, 
have been proposed (see, for example, refs 0], [H] and 0]). The theoretical 
predictions are however in disagreement with the experiments and this is 
mainly due to the fact that, in the series expansion, the terms of superior 
order are greater than those of inferior order. Therefore truncation of the 
series at some order is not mathematically justified. 

A thermodynamic field theory (TFT) has been proposed in 1999 in order to 
evaluate how the relations between fluxes and forces, Eqs (pQ), deform when 
the thermodynamic system is far from the linear (" Onsager" ) region [7|] . At- 
tempts to derive a generally covariant thermodynamic field theory (GTFT) 
can be found in refs 0|. The characteristic feature of the TFT is its purely 
macroscopic nature. This does not mean a formulation based on the macro- 
scopic evolution equations, but rather a purely thermodynamic formulation 
starting solely from the entropy production and from the transport equa- 
tions, i.e., the fluxes-forces relations. The latter provide the possibility of 
defining an abstract space (the thermodynamic space), covered by the n 
independent thermodynamic forces X^, whose metric is identified with the 
symmetric part of the transport matrix. The law of evolution is not the dy- 
namical law of particle motion, or the set of two-fluid macroscopic equations 
of plasma dynamics. The evolution in the thermodynamic space is rather 
determined by postulating three purely geometrical principles: the Shortest 
Path Principle, the Thermodynamic Field Strength in closed form, and the 
Principle of Least Action. From theses principles, a set of field equations, 
constraints, and boundary conditions are derived. These equations, referred 
to as the thermodynamic field equations, determine the nonlinear corrections 
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to the linear ("Onsager") transport coefficients. However, the formulation of 
the thermodynamic field theory, as reported in refs 0], raises the following 
fundamental objection 

There are no strong experimental evidences supporting the requirement that 
the thermodynamic field strength is in a closed form. 

Moreover, the principle of general covariance, which in refs [8[ has been as- 
sumed to be valid for general transformations in the space of thermodynamic 
configurations, is, in reality, respected only by a very limited class of thermo- 
dynamic processes. In this paper, through an appropriate mathematical for- 
malism, the iso-entropic formalism, the entire TFT is re-formulated remov- 
ing the assumptions regarding the closed-form of the thermodynamic field 
strength and the general covariance principle (GCP). The GCP is replaced 
by the thermodynamic covariance principle (TCP), or the Prigogine state- 
ment establishing that thermodynamic systems, obtained by a trans- 
formation of forces and fluxes in such a way that the entropy production 
remains unaltered, are thermodynamically equivalent. This principle ap- 
plies to transformations in the thermodynamic space, and they are referred 
to as the thermodynamic coordinate transformations (TCT). It is worth- 
while mentioning that the TCP is actually largely used in a wide variety of 
thermodynamic processes ranging from non equilibrium chemical reactions 
to transport processes in tokamak plasmas (see, for examples, the papers 
and books cited in refs [1]) and [13]). To the author knowledge, the valid- 
ity of the thermodynamic covariance principle has been verified empirically 
without exception in physics until now. 

The analysis starts from the following observation. Consider a relaxation 
process of a thermodynamic system in the Onsager region. If the system 
relaxes towards a steady-state along the shortest path in the thermodynamic 
space, then the Universal Criterion of Evolution is automatically satisfied. 
Indeed, in this case, we can write 



where the dot over the variables indicates the derivative with respect to the 
arc parameter <j, defined as 



Parameter q can be chosen in such a way that it vanishes when the system 
begins to evolve and it assumes the value, say I, when the system reaches 
the steady-state. In the Onsager region, the thermodynamic space is an 



(6) 



df = (L^dX»dX») 1/2 



(7) 



4 



Euclidean space with metric L^ u . The equation of the shortest path reads 
X^ = 0, with solution of the form 

X» = + V 1 (8) 

where and 6 M are arbitrary constant independent of the arc parame- 
ter. Inserting Eq. (jSJ) into Eq. © and observing that L^a^a? = 1 and 
fo^a^a" = 0, we find 

J^X* = q + To^aW (9) 

At the steady state (i.e. for q = I) J^X^ \ st . s tate= (because V \ st . s tate= 0). 
Eq. ([9]) can then be written as 

p = -(I - q ) < (with P = J^) (10) 



or 

"—('-0(^-5-5-) S« (H) 

The equation for the dissipative quantity P, when the thermodynamic sys- 
tem relaxes in the linear region, is thus given by Eq. 0: 

Also note that & = 2P < i.e., the minimum entropy production theorem 
is also satisfied during relaxation. Outside of Onsager's region, one may be 
tempted to construct a Riemannian space (of 3 or more dimensions) which is 
projectively flat i.e., having a vanishing Weyl's projective curvature tensor. 
In this case, indeed, there exists a coordinate system such that the equations 
of the shortest path are linear in the coordinates [i.e., the shortest paths are 
given by equations of the form (jSJ)]. In this respect, we have the following 
Weyl theorem [lj]]: a necessary and sufficient condition that a Riemannian 
space be projectively flat is that its Riemannian curvature be constant ev- 
erywhere. On the other hand, to re-obtain the Onsager relations, we should 
also require that, near equilibrium, the Riemannian space reduces to a flat 
space (which has zero Riemannian curvature). The Weyl theorem can be 
conciliated with our request only if there exists a coordinate system such 
that Eqs (|2|) are valid everywhere, which is in contrast with experiments. 
Thus one wants the Universal Criterion of Evolution satisfied also out of 
the Onsager region, without imposing a priori any restrictions on transport 
coefficients, a non- Riemannian thermodynamic space is required. Clearly, a 
transport theory without a knowledge of microscopic dynamical laws can not 
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be developed. Transport theory is only but an aspect of non-equilibrium sta- 
tistical mechanics, which provides the link between micro and macro-levels. 
This link appears indirectly in the "unperturbed" matrices, i.e. the L^ v and 
the fo^u coefficients, used as an input in the equations: these coefficients 
have to be calculated in the usual way by kinetic theory. 
In section [21 we introduce a non-Riemannian space whose geometry is con- 
structed in such a way that 

A. The theorems valid when a generic thermodynamic system relaxes out 

of equilibrium are satisfied] 

B. The differential equations for the transport coefficients are covariant 
under the thermodynamic coordinate transformations (TCT). 

We shall see that the properties of geometry do not depend on the short- 
est paths but upon a particular expression of the affine connection. Our 
geometry is then of affine type and not of projective type. At the end of 
section [2] we derive the field equations for the transport coefficients through 
an appropriate mathematical formalism: the iso-entropic formalism. This 
formalism allows to respect the Prigogine statement. New objects like ther- 
modynamic covariant differentiation or the thermodynamic curvature are 
also introduced. We shall see that under the weak-field approximation and 
when <7 3> 1, but only in these limits, the new thermodynamic field equa- 
tions reduce to the ones obtained in refs [7(. So that, all results found in 



refs [12(|, for magnetically confined plasmas, and in refs 13], for the nonlin- 
ear thermoelectric effect and the unimolecular triangular reaction, remain 
valid. In section [3] it is shown that this formalism is able to verify the ther- 
modynamic theorems (in particular, the Universal Criterion of Evolution) 
for systems relaxing out of the Onsager region. Mathematical details and 
demonstrations of the theorems are reported in the annexes. 

2 The Iso-Entropic Formalism 

Consider a thermodynamic system driven out from equilibrium by a set of n 
independent thermodynamic forces {X^} (/U = 1, • • • n). It is also assumed 
that the system is submitted to time-independent boundary conditions. The 
set of conjugate flows, { J^}, is coupled to the thermodynamic forces through 
the relation 

J M = r„ u X v (13) 
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where denote the transport coefficients. The symmetric piece of t^ v is 
denoted with g^ u and the skew-symmetric piece as fav- 

where 

ft/At ( 15 ) 

-/^ (16) 

It is assumed that g^i, is a positive definite matrix. With the elements 
of the transport coefficients two objects are constructed: operators, which 
may act on thermodynamic tensorial objects and thermodynamic tensorial 
objects, which under coordinate (forces) transformations, obey to well spec- 
ified transformation rules. 

Operators 

Two operators are introduced, the entropy production operator a {X) and 
the dissipative quantity operator P(X), acting on the thermodynamic forces 
in the following manner 

a(X) :-» <r(X) = XgX T 

P(X) :-► P(X) = XtX t (17) 

In Eqs (|17p . the transport coefficients are then considered as elements of 
the two n x n matrices, r and g. The positive definiteness of the matrix 
g^y ensures the validity of the second principle of thermodynamics: a > 0. 
These matrices multiply the thermodynamic forces X expressed as n x 1 
column matrices. The dot symbol stands for derivative with respect to 
parameter defined in Eq. (|29l) . Thermodynamic states X s such that 

P(X S ) = (18) 

are referred to as steady-states. These are physical quantities and should re- 
main invariant under thermodynamic coordinate transformations. Eqs (|17p 
should not be interpreted as the metric tensor g^ u , which acts on the coor- 
dinates. The metric tensor acts only on elements of the tangent space (like 
dX^, see the forthcoming paragraphs) or on the thermodynamic tensorial 
objects. 

Transformation Rules of Entropy Production, Forces and Flows 



I flU 



- IT, 



fj,u 



Tun) 
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According to Prigogine's statement [lj, thermodynamic systems are ther- 
modynamically equivalent if, under transformation of fluxes and forces the 
bilinear form of the entropy production, a, remains unaltered. In mathemat- 
ical terms, this implies: 

a = J a X» = J'^X"* (19) 

This condition requires that the transformed thermodynamic forces and 
flows satisfy the relation 

X* = ^—X v 
dX v 

pi vv 

J '» = wrS* ( 2 °) 

These transformations are referred to as Thermodynamic Coordinate Trans- 
formations (TCT). The expression of entropy production becomes accord- 
ingly 

a = J^X* = t^X^X" = 9liV X»X v = g'^X'^X'" = a' (21) 
From Eq. (|2ip we find 

, __ dX» dX u 
9\ K - V^pl^ ( 22 ) 

Moreover, inserting Eqs ([20]) and Eq. ([22]) into relation J M = (g^ v + f^X", 
we obtain 

r-( a > +f dx » dxv \ x *K m . 



dX' x dX 
or 

J'x = (gL + fL)x' K with fL = U» dx , Xdx/K (24) 

Hence, the transport coefficients transform like a thermodynamic tensor of 
second orderQ 

Properties of the TCT 

By direct inspection, it is easy to verify that the general solutions of equa- 
tions ([20]) are 

^^(fr.^-^) (*> 



2 We may qualify as thermodynamic tensor (taken as a single noun) a set of quantities 
where only transformations Eqs \2U\) are involved. This is in order to qualify as a tensor, 
a set of quantities, which satisfies certain laws of transformation when the coordinates 
undergo a general transformation. Consequently every tensor is a thermodynamic tensor 
but the converse is not true. 
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where is an arbitrary function of variables X J /X^ with (j = 2, . . . , n). 
We may (or we must) require that, in the Onsager region, Eqs (|25p reduce 
to linear homogeneous transformations 

X'» = ($X U (26) 

where Cy are constant coefficients (i.e., independent of the thermodynamic 
forces). Note that from Eq. (l20j) . the following important identities are de- 
rived 

d 2 X'v d 2 X fi 
yy u _ — n • y' v ——— — n (r>7\ 

A dx-dx-~ u ' dx'-dx'-~ u ( <) 



Moreover 



d dX v d 



dX'» dX'v dX v ^ 

i.e., dX^ and d/dX^ transform like a thermodynamic contra- variant and 
a thermodynamic covariant vector, respectively. According to Eq. (|28p. 
thermodynamic vectors dX^ define the tangent space to Ts. It also follows 
that the operator P(X), i.e. the dissipation quantity, and in particular the 
definition of steady-states, are invariant under TCT. Parameter defined 
as 

df = g^dXW (29) 
is a scalar under TCT. The operator O 

O = X»—- = X'^—— = O' (30) 
dX^ dX'v K ' 

is also invariant under TCT. This operator plays an important role in the 
formalism. 

Thermodynamic Space, Thermodynamic Covariant Derivatives and 
Thermodynamic Curvature 

A non-Riemannian space with a linear connection T^g is now introduced. 
Consider an n-space in which the set of quantities T^ g is assigned as func- 
tions of the n independent thermodynamic forces X^, chosen as coordinate 
system. Under a coordinate (forces) transformation, it is required that the 
functions transform according to the law 



dX'v dX x dX K dX'» d 2 X u 
ap Xk dX» dX«* dX'P dX" dX la dX'P 
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With the linear connection r£°g, the absolute derivative of a thermodynamic 
contra-variant vector along a curve can be defined as 

T = T + r "» r T (32) 

It is easily checked that, if the parameter along the curve is changed from 
<j to g, then the absolute derivative of a thermodynamic tensor field with 
respect to g is ck/dg times the absolute derivative with respect to The 
absolute derivative of any contra- variant thermodynamic tensor may be eas- 
ily obtained generalizing Eq. (I32p . In addition, the linear connection is 
submitted to the following basic postulates: 

1. The absolute derivative of a thermodynamic contra-variant tensor is a 

thermodynamic tensor of the same order and type. 

2. The absolute derivative of an outer product of thermodynamic tensors, 

is given, in terms of factors, by the usual rule for differentiating a 
product. 

3. The absolute derivative of the sum of thermodynamic tensors of the 
same type is equal to the sum of the absolute derivatives of the ther- 
modynamic tensors. 

In a space with a linear connection, we can introduce the notion of the 
shortest path defined as a curve such that a thermodynamic vector, initially 
tangent to the curve and propagated parallelly along it, remains tangent to 
the curve at all points. By a suitable choice of the parameter g, the differ- 
ential equation for the shortest path is simplified reducing to 

d 2 X^ dX a dX? _ 

~lg^ +T ^^g~^g~~° (33) 

To respect the general requirement A. (see section [3]), it is required that the 
absolute derivative of the entropy production satisfies the equality 

More in general, it is required that the operations of contraction and absolute 
differentiation commute for all thermodynamic vectors. As a consequence, 
the considered space should be a space with a single connection. The ab- 
solute derivative of a covariant thermodynamic vector T^ is then defined 
as 

5T U dT u „ dX? , N 

'st-^-sr (35) 
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The absolute derivative of the most general contra- variant, covariant and 
mixed thermodynamic tensors may be obtained generalizing Eqs (j32j) and 

The derivatives, covariant under TCT, of thermodynamic vectors, are de- 
fined as 

W dX v av 
dT 

For the entropy production, it is also required that 

P^v = 0>|(U ( 37 ) 



More in general, Eq. (|37p should be verified for any thermodynamic scalar 

a/3 



T. This postulate requires that the linear single connection T^r, is also 



symmetric i.e., T^g = T^ a . A non-Riemannian geometry can now be con- 
structed out of n 2 (n + l)/2 quantities, the components of T^, according to 
the general requirements A and B mentioned in the introduction. 
In the forthcoming paragraph, the expression of the affine connection T^a is 
determined from assumption A. In section [3] it is shown that the Universal 
Criterion of Evolution, applied to thermodynamic systems relaxing towards 
a steady-state, is automatically satisfied along the shortest path if, in case 
of symmetric processes (i.e., for a — a processes), we impose 

r ^ - 2 g W + dx^~ dx*) + Ya X °^ (38) 

In the general case, we have 
T ap - N» g KX | a/? | + —X K 0{g aP ) + —X K X + — ) 

+^-f Ki X X + —) (39) 

where the thermodynamic Christoffel symbols of the second kind are intro- 
duced 

apj 2 g \dX?^dX a dX x ) V ' 

and matrix N^ K is defined as 

Ny* = 9»v + ^f„ K X K X v + ^f VK X K X„ with : N^N UK = 8$ (41) 
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In appendix [6] it is proven that the affine connections Eqs (|38|) and (|39p 
transform, under a TCT, as in Eq. (|3ip and satisfy the postulates 1., 2. and 
3. From Eq. (|4ip we obtain 



N^ V X U = yg^u + —f^nX^Xy + —f VK X K X^X y = 

N^X" = [g^ + ^ + = J, (42) 

N^X* = = a 

While 

jv^% = N^N^X V = N K »N Vfl X u = X K (43) 
N^ K J K = N K »J K = X^ 
N^J K J, L =X"J li = a 

The shortest path is the same for two symmetric connections whose coeffi- 
cients are related as 

f^ = r^ + ^/j + ^« (44) 

where ip a is an arbitrary covariant thermodynamic vector and by denotes the 
Kronecker tensor. In literature, modifications of the connection similar to 
Eqs (|44p are referred to as projective transformations of the connection and 
ip a the projective covariant vector. The introduction of the affine connection 
gives rise to the following difficulty: the Universal Criterion of Evolution is 
satisfied for every shortest path constructed with affine connections T^o, 
linked to T^a by projective transformations. This leads to an indetermi- 
nation of the expression for the affine connection, which is not possible to 
remove by using the Prigogine statement and the thermodynamic theorems 
alone. This problem can be solved by postulating that the thermodynamic 
field equations (i.e., the field equations for the affine connection and the 
transport coefficients) be symmetric and projective-invariant (i.e., invariant 
under projective transformations). 

For any covariant thermodynamic vector field T u , we can form the thermo- 



dynamic tensor R^ Xk in the following manner [14| 



T v \x\ K ~ T„|«|A = t hK*k ( 45 ) 

where 

fir'* r?r M 

pj U _ UL UK __ W± y\ j_ r l) r /l _ p77 -pfj, ( A o\ 

^uXk — Q X \ Q-gK VK r)X 1 vX L r]K V* 1 ) 
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with R^\ K satisfying the following identities 



_ _ dm 

Kxn + Kku + R Xun = ( 47 ) 

#\ , =0 

VAK\r\ VKT]\X Vf]X\K 

By contraction, we obtain two distinct thermodynamic tensors of second 
order 

**> - 2 V- 2 UlA Q X u) ^) 

with F\ v being skew-symmetric and R u \ asymmetric. Tensor R u \ can be 
re-written as 

RuX = B vX + Fxu where 



-B^A - #Ak — 2 ^ A + J QX^ + r ^ r »? A r ^A r w ( 49 ) 

Hence, i^x^ is the skew-symmetric part of i?^A- It is argued that the thermo- 
dynamic field equations can be derived by variation of a stationary action, 
which involves R v \- Symmetric and projective- invariant field equations may 
be obtained by proceeding in following manner: 1) a suitable projective 
transformation of the affine connection is derived so that R u \ be symmetric 
and F\ v be a zero thermodynamic tensor and 2) the most general projective 
transformation that leave unaltered R u \ and F\ v (= 0) is determined. By a 
projective transformation, it is found that 



B v x = B vX + n 



^ = ^ + -2-^-^ (50) 

Eq . (BSD shows that F A; , can be written as the curl of the vector a„/2 defined 
as [l5|] 

^ = 1^-1*1 (51) 



Consequently, by choosing 



*» = ^rfc-LlM (52) 



n + 1 V I 
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we have F\ v = and R u \ = B v \ . From Eqs ([50]) , we also have that the ther- 
modynamic tensor R v \ remains symmetric for projective transformations of 
connection if, and only if, the projective covariant vector is the gradient of 
an arbitrary function of the X's [15J. In this case, the thermodynamic tensor 
F\ v remains unaltered i.e., F\ v = 0. Hence, at this stage, the expression 
of the affine connection is determined up to the gradient of a function, say 
6, of the thermodynamic forces, which is also scalar under TCT. Let us im- 
pose now the projective-invariance. Eqs (|50p indicate that a necessary and 
sufficient condition that R u \ be projective-invariant is that 

d 2 d> 86 86 



8X x dX» dX x 8X» 
6 = 



with 



d 2 <f> 
dX^dX v 



(in the Onsager region) (53) 





where 6 is a function, invariant under TCT. The solution of Eq. (|53|) is 6 = 
everywhere. The final expression of the affine connection for symmetric 
processes reads then 

The general case is given by 



8X^ 



where 

, N^g KX f A } N""X K N"« vX fdf v x , Of, 

n + 1 \w) 2(n+l)a {9vr,) 2(n + l)a K \dX u dX 

2(n + l)a U " X X \8X^ + 8W ) + n + 1 ^ (5bj 

Note that the thermodynamic space tends to reduce to a (thermodynamic) 
Riemannian space when a v -C 1. The following definitions are adopted: 

• The space, covered by n independent thermodynamic forces X^, with 
metric tensor and a linear single connection given by Eq. (|55|) , may 
be referred to as thermodynamic space Ts (or, space of the thermody- 
namic forces). 
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In Ts, the length of an arc is defined by the formula 



L = 




dX^ dX v \ 1/2 




(57) 



The positive definiteness of matrix g^ v ensures that L > 0. Consider a 
coordinate system X^, defining the thermodynamic space Ts. 

• All thermodynamic spaces obtained from Ts by a TCT transformation, 
may be called iso-entropic spaces. 

In the TFT description, a thermodynamic configuration corresponds to a 
point in the thermodynamic space Ts. The equilibrium state is the origin of 
the axes. Consider a thermodynamic system out of equilibrium, represented 
by a certain point, say a, in the thermodynamic space 

• A thermodynamic system is said to relax towards another point of 
the thermodynamic space, say b, if it moves from point a to point b 
following the shortest path (133f) . with the affine connection given in 



• With Eq. (|55p . Eqs ()36[) may be called the thermodynamic covariant 
differentiation of a thermodynamic vector while Eqs. (I32|) and (I35p the 
thermodynamic covariant differentiation along a curve of a thermody- 
namic vector. 

• With affine connection Eq. (|55p . R^\ K may be called the thermody- 
namic curvature tensor. 

• The scalar R obtained by contracting the thermodynamic tensor R v \ 
with g uX (i.e. R = R u \g uX ) may be called the thermodynamic curva- 
ture scalar. 

The Principle of Least Action 

From expression (|55p , the following mixed thermodynamic tensor of third 
order can be constructed 



Eq. d55D- 



a/3 — 



{ 



A 

a(3 




2a 



A 



(§n+ft)+^ + « -{:,} < 58 > 



15 



This thermodynamic tensor satisfies the important identities 

Kp = Kp = ° ( 59 ) 

Again, from the mixed thermodynamic tensor of fifth order can be 
constructed 

sZ P - \ (Va + Kx% + + Kx% - - Kpti) (so) 

By contraction, a thermodynamic tensor of third order, a thermodynamic 
vector and a thermodynamic scalar can be formed as follows 

(61) 

The following postulate is now introduced: 

There exists a thermodynamic action I, scalar under TCT , which is sta- 
tionary with respect to arbitrary variations in the transport coefficients and 
the affine connection. 

This action, scalar under TCT, constructed from the transport coefficients 
and their first and second derivatives, should have linear second derivatives. 
In addition, it should be stationary when the affine connection takes the 
expression given in Eq. (|55j> . The only action satisfying these requirements 
is 

gTy/gdTX (62) 

where d n X denotes an infinitesimal volume element in Ts and T*^ is the 
expression given in Eq. ([55]) i.e., = ty* v + i/JlV To avoid misun- 
derstanding, while it is correct to mention that this postulate affirms the 
possibility of deriving the thermodynamic field equations by a variational 
principle it does not state that the expressions and theorems obtained from 
the solutions of the thermodynamic field equations can also be derived by 
a variational principle. In particular the well-known Universal Criterion 
of Evolution established by Glansdorff-Prigogine can not be derived by a 
variational principle (see also section [3]). 

The Thermodynamic Field Equations 



R 
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As a first step, the transport coefficients and the affine connection are 
subjected to infinitesimal variations i.e., — ► g^ + 5g^ u , f^ v — ► + 
and -> + tfT^, where Sg^, Sf^ and are arbitrary, 

except that they are required to vanish as | |— > oo. By imposing that 
the action (|62p is stationary with respect to arbitrary variations in g^, 
fnv and r* the thermodynamic field equations (i.e., the equations for the 
transport coefficients and the affine connection) are derived. The results 
are, respectively (see appendix [7]) 

i r)f x 
~~ ~~5g^ 

Sf — ^ = (63) 

where the variations of the affine connection (|55p with respect to the trans- 
port coefficients appear in the first two equations. From the first equation 
of Eqs (I63h , the expression for the thermodynamic curvature scalar is given 
as 

9 sr A 

R = —g^sf-^- (64) 



The third equation of Eqs (|63p can be re-written as 

5^,A - ^\9au ~ Fv\9a» = ~^\9au ~ ^u\9ocn (65) 

where the comma (, ) denotes partial differentiation. Adding to this equation 
the same equation with \x and A interchanged and subtracting the same 
equation with v and A interchanged gives 

9»v,\ + 9\v,n ~ 9fi\,u = ^9av^xa - Igav^^ (66) 



or 



n^={^}+n^=K (67) 

For a — a processes, close to the Onsager region, it holds that 
g»v = + h^v + 0(e 2 ) 

X a = 0(e) with e = Ma J \ E ^nvalues[g -L ]\ , <i 

I Eigenvalues[L „„J J 
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where A CT = 1/a and h^ u are small variations with respect to Onsager's 
coefficients. In this region, Eq. (|62p is stationary for arbitrary variations of 
h/jj, and r^,. It can be shown that [7|] 

tXk d 2 /y xk d<1}l \n _ T x K d 2 h Xv _ Xk d 2 h Xfl = 2 
dX x dX* dX^dX" dX^dX^ dX^dX" K ' 

9^,x = + O(e 2 ) or Y% = h^(h m<u + - h^ v ) + 0(e 2 ) (69) 

where the semicolon indicates the covariant derivative performed with the 
Christoffel symbols. Eqs (i68j) should be solved with the appropriate gauge- 
choice and boundary conditions. 

The validity of Eqs (I69p has been largely tested by analyzing several sym- 
metric processes, such as the thermoelectric effect and the unimolecular 
triangular chemical reactions [7|. More recently, these equations have been 
also used to study transport processes in magnetically confined plasmas. In 
all examined examples, the theoretical results of the TFT are in line with 
experiments. It is worthwhile mentioning that, for transport processes in 
tokamak plasmas, the predictions of the TFT for radial energy and matter 
fluxes are much closer to the experimental data than the neoclassical the- 
ory (based on linearized Boltzmann's equation), which fails with a factor 
10 3 -=- 10 4 [u} and 0. 

Some Remarks on Spatially Extended Thermodynamic Systems 

The macroscopic description of thermodynamic systems gives rise to 
state variables that depend continuously on space coordinates. In this case, 
the thermodynamic forces possess an infinity associated to each point of the 
space coordinates. The system may be subdivided into N cells (NxNxN in 
three dimensions), each of which labeled by a wave-number k, and we follow 
their relaxation. Without loss of generality, we consider a thermodynamic 
system confined in a rectangular box with sizes l x , l y and l z . Denoting with 
k the wave-number 

n x = 0,±1,--- ±N X 
2tt(^,^,^) with { n y = 0,±l,---±N y (70) 



111 

L z 



n z = 0,±l,---±N z 
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the fluxes and forces, developed in (spatial) Fourier's series, read 

N 



n=-N 
N 



X»(r,t)=J2 lf k0 (i)exp(ik'.r) (71) 



L (k') 

n'=-N 



where, for brevity, n and N stand for n = (n x ,n y , n z ) and N = (N x , N y , N z ), 
respectively. The Fourier coefficients are given by 

^u(k) (*) = 7^ J Jv( r , t) exp(-ik • r)dv 

Xfa® = ^J^(r,t)e W (-fc'-r)dv (72) 

In particular, the contributions at the thermodynamic limit (i.e., for k — > 0) 
are expressed as 

X{ 0) (t) = ~J X> t (T,t)dv=X>>(t) (73) 

The entropy production and the transport relations take, respectively, the 
form 

a(r,t) = J tl (r,t)X»{r,t) > 

Mr,t)=r^(r,t)X»(r,t) (74) 

Considering that 

fix fly flz 

/ / / exp[i(k + k') • r]dv = 0, <5k+k',o with (75) 
Jo Jo Jo 



, if k + k' + _ 777 

Ok+k',o = \ ^ ifk + k' — a = y 



from the first equation of Eq. (|74p we also find 

f J,(r,t)X»(r,t) ^ = o(j M0) (t)lf 0) (i) + ^J Kk) (t)lf_ k) (t)) >0 



(76) 
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On the other hand, we have 

Jm (*) = V( ) (*)*(o) (*) + E V(k) (*)*(-k) (*) (77) 

k^O 

where 

V(k)(*) = ^ J r flu (r,t)exp(-ik ■ r)dv (78) 
Eq. (|76p can then be brought into the form 



a dv = Qg^(t)Xfo(t)Xfa(t) 



k^O 

where 

9^(k)(t) = J q G^( r ' *) exp(-ik • r)tfu with 

= ^V^M) + T ^M)] (80) 

In Eq. (|79p . the first term is the contribution at the thermodynamic limit 
whereas the second expression reflects the interactions between the k-cell 
and the other cells. In a relaxation process, contributions from different 
wave-numbers are negligible with respect to those with same wave-numbers 
and, hence, we finally obtain 

f <rdv~ fi^ K0) (t)Xg } (t)Xfo (t)>0 V XJ 3 (i) (and a # 0) (81) 
Last inequality is satisfied /or any XiL (i) if, and only if 

9^u(o){t) = Q^ir.t) dv = g^ v {t) (82) 

is a positive definite matrix. Therefore, for spatially extended thermody- 
namic systems, we replace X^it) — ► -X^(i) and T^ v {t) — > r^^it). Un- 
der these conditions, Eqs (|63p determine the nonlinear corrections to the 
Onsager coefficients while Eqs (j32|) and Eqs (|36p. with affine connection 
Eq. (|55p . are the thermodynamic covariant differentiation along a curve and 
the thermodynamic covariant differentiation of a thermodynamic vector, re- 
spectively. 
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The Privileged Thermodynamic Coordinate System 

By definition, a thermodynamic coordinate system is a set of coordinates 
defined so that the expression of the entropy production takes the form of 
Eq. (|19[) . Once a particular set of thermodynamic coordinates is determined, 
the other sets of coordinates are linked to the first one through a TCT [see 
Eqs (|20p ]. The simplest way to determine a particular set a coordinate 
transformation is to quote the entropy balance equation 

^+V-3 s = a (83) 
at 

where ps is the local total entropy per unit volume and J s is the entropy 
flux. Let us consider, as an example, a thermodynamic system confined in a 
rectangular box where chemical reactions, diffusion of matter, macroscopic 
motion of the volume element (convection) and heat current take place si- 



multaneously. The entropy flux and the entropy production read [la] , [171] 

J s = — (J q - ^ Jj/Xj) + y^pjVjSi 



T 

- = -V^4E [V(f H +£^4£ n *A^ > (84) 



(85) 



where pi, piS{ and A{ are the chemical potential, the local entropy and the 
affinity of species "i", respectively. Moreover, Fj indicates the external force 
per unit mass acting on Iljj the components of the dissipative part of 
the pressure tensor Aiij (M-ij = p5ij + IL^-; p is the hydrostatic pressure) 
and Vj is the component of the hydrodynamic velocity (see, for example, ref. 
[l~8| ) . The set of thermodynamic coordinates is given as 



'v-- -- 

> T' ~T 



For this particular example, this set may be referred to as the privileged 
thermodynamic coordinate system. Other examples of privileged thermody- 
namic coordinate system, concerning magnetically confined plasmas, can be 
found in refs Q, O and 0. 
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3 Thermodynamic Theorems for Systems out of 
Equilibrium 



In 1945-1947, Prigogine proved the minimum entropy production theorem 
[l|, which concerns the relaxation of thermodynamic systems near equilib- 
rium. This theorem states that: 

Minimum Entropy Production Theorem (MEPT) 

For a — a processes, a thermodynamic system, near equilibrium, relaxes 
to a steady-state X s in such a way that the inequality 

% <- ° < 87 > 

is satisfied throughout the evolution and is only saturated at X s . 
The minimum entropy production theorem is generally not satisfied far from 
equilibrium. However, P. Glansdorff and I. Prigogine demonstrated in 1964 
that a similar theorem continues to hold for any relaxation to a steady-state, 
which reads 0] 

Universal Criterion of Evolution (UCE) 

When the thermodynamic forces and conjugate flows are related by a 
generic asymmetric tensor, regardless of the type of processes, for time- 
independent boundary conditions a thermodynamic system, even in strong 
non- equilibrium conditions, relaxes towards a steady-state in such a way that 
the following universal criterion of evolution is satisfied: 

dX^ 

V S ./„— < (88) 
This inequality is only saturated at X s . 

For a — a processes, the UCE reduces to the MEPT in the Onsager region. 
Again, Glansdorff and Prigogine demonstrated this theorem using a purely 
thermodynamical approach. In this section we shall see that if the system 
relaxes towards a steady-state along the shortest path then the Universal 
Criterion of Evolution is automatically satisfied. 

By definition, a necessary and sufficient condition for a curve to be the 
shortest path is that it satisfies the differential equation 

d 2 X» „ dX a dX? . . dX* 

where (p(t) is a determined function of time. If we define a parameter g by 
dQ f ... » nl dX» 



dt 



f * , dX v . 

cexp ip dt with ip = ip — 2tp u (90) 
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where c is an arbitrary constant and ip v the projective covariant vector, 
Eq. ([89]) reduces to Eq. ([33]) with T^a given by Eq. ([39]) . Parameter g is not 
the affine parameter s of the shortest path. The relation between these two 
parameters is 

g = b j exp(-2 J ip u dX u ^Jds (91) 



where b is an arbitrary constant. Eq. (|90l) allows us to choose the parame- 
ter g in such a way that it increases monotonically as the thermodynamic 
system evolves in time. In this case, c is a positive constant and, without 
loss of generality, we can set c = 1. Parameter g can also be chosen so 
that it vanishes when the thermodynamic system begins to evolve and it 
takes the (positive) value, say I, when the system reaches the steady-state. 
Multiplying Eq. (|33p with the flows and contracting, we obtain 

d 2 X^ u dX a dX 13 , . 

h-inr + JA , =0 (92) 



dg 2 af3 dg dg 



However 



d 2 X» = dP_ /*\ 2 _ dX a dX? x 8g aX _ dX a dX? x 8f aX 



M dg 2 ~ dg \dg) dg dg dX? dg dg dX? K ' 

where P = J^^jf an d after taking into account the identities fixv^j—^f- 
( ) and = 1- I n addition, recalling Eq. ([43 p and the relations 

X^X^ = a and f flu X tJ, X u = 0, it can be shown that 

dX^dX^ = dX^dX^ x dg aX dX a dX\ x df aX 
«P dg dg dg dg OX? + dg dg OX? [ ' 



Summing Eq. f)93f) with Eq. (J94J) and considering Eq. (|92p . gives 

(95) 

Integrating Eq. ([9"5]) ] from the initial condition to the steady-state, we find 



dP _ /d?\2 
dg \dgJ 



PCVJ-P /(^i) 2 ^>0 (96) 



^dgJ 

From Eq. (j!8p we have P(X S ) = P{X s )dq/dg = 0, so we finally obtain 



P = J ft — = - / -1 d^<0 (97) 
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where the inequality is only saturated at the steady-state. Recalling Eq. ([90]) . 
the inequality established by the UCE can be derived 

V = pd J t = J ^ ( ex P / ( 98 ) 

Eq. ([95]) can be re-written as 



d_ 



dg) 



P 



dg 



(99) 



This equation generalizes Eq. (|12p . which was valid only in the near equi- 
librium region (notice that, in the linear region, dq/dg = 1/6). Integrating 
Eq. (|99p . the expression of the dissipative quantity P is derived 

\dq) _/ f V dg ) v ^ dg dg ) J q V ^ dg dg J ^ 

(100) 

On the right, it is understood that the X's are expressed in terms of (){$)■ 
Eq. (jlOOp generalizes Eq. (|10p . which was valid only in the linear region. For 
a — a processes in the Onsager region, Eq. (|98p implies the validity of the 
inequality (fBTl) . Indeed, Eq. (f34l) gives 

" = J ^ + ^ = 2J ^ + ^A (101) 

d£ dg dg dg 

In the linear region, the coefficients of the affine connection vanish. Eq. ()101|) 
is simplified reducing to 



da da dg / dX^ 1 c?£> 
(it dg dt V ^ <ig 



2( J^-j— 37) = IV < (102) 



where the inequality is saturated only at the steady state. Let us now 
consider spatially extended thermodynamic systems. In space-dependent 
systems, the dissipative quantity should be expressed in the integral form 

V= [ J^r,t)dtX^(r,t) dv (103) 
Jn 

where diX^ = dX 11 /dt. In terms of wave- vectors k, Eq. (|103p can be brought 
into the form 

v = n(j m (t)dtxt 0) (t) + ^ J^m^-v®) (io4) 

k^0 
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where Eq. (|75p has been taken into account. As already mentioned in section 
[2j in a relaxation process, contributions from different wave-numbers are 
negligible with respect to those with same wave-numbers and, hence, we 
finally obtain 

V ~ n,J^{t)d t X^{t) < (105) 

where Eq. (|98p has been taken into account. It is therefore proven that 
the Universal Criterion of Evolution is automatically satisfied if the system 
relaxes along the shortest path. Indeed it would be more exact to say: the 
affine connection, given in Eq. (j39|) . has been constructed in such a way 
that the UCE is satisfied without imposing any restrictions on the trans- 
port coefficients (i.e., on matrices g^ v and f^ u ). In addition, analogously to 
Christoffel's symbols, the elements of the affine connection are constructed 
from matrices g^ v and and their first derivatives in such a way that all 
coefficients vanish in the Onsager region. Eq. (|39|) provides the simplest 
expression satisfying these requirements. 

The Minimum Rate of Dissipation Principle (MRDP) 

In ref.[s| the validity of the following theorem is shown: 
The generally covariant part of the Glansdorff-Prigogine quantity is always 
negative and is locally minimized when the evolution of a system traces out 
a geodesic in the space of thermodynamic configurations. 
It is important to stress that this theorem does not refer to the Glansdorff- 
Prigogine expression reported in Eq. ([88]) but only to its generally covariant 
part. Moreover, it concerns the evolution of a system in the space of ther- 
modynamic configurations and not in the thermodynamic space. One could 
consider the possibility that the shortest path in the thermodynamic space 
is an extremal for the functional 

J„X^ (106) 



SI 



The answer is negative. Indeed, a curve is an extremal for functional 
Eq. (|106p if, and only if, it satisfies Euler's equations^] 

As it can be easily checked, this extremal coincides with the shortest path 
if 

2\dXP dX a , 



2{dX^ + dX^)- T ^ M ^ = ° Wh6re (1 ° 8) 
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and r^g given in Eq. ([55]) , However, Eqs. (|108j) are n 2 (n + l)/2 equations 
for n 2 variables (the transport coefficients) and, in general, for they 
do not admit solutions. We have thus another proof that the Universal 
Criterion of Evolution can not be derived from a variational principle. 

4 Conclusions and Limit of Validity of the Ap- 
proach 

The main purpose of this paper is to present a new formulation of the 
thermodynamic field theory (TFT) where one of the basic restrictions in the 
old theory, the closed- form of the thermodynamic field strength (see ref . [?J ) , 
has been removed. Furthermore, the general covariance principle, respected, 
in reality, only by a very limited class of thermodynamic processes, has 
been replaced by the thermodynamic covariance principle, first introduced 
by Prigogine for treating non equilibrium chemical reactions. The validity 
of the Prigogine statement has been successfully tested, without exception 
until now, in a wide variety of physical processes going beyond the domain of 
chemical reactions. The introduction of this principle requested, however, 
the application of an appropriate mathematical formalism, which may be 
referred to as the iso-entropic formalism. The construction of the present 
theory rests on two assumptions: 

• The thermodynamic theorems valid when a generic thermodynamic 
system relaxes out of equilibrium are satisfied; 

• There exists a thermodynamic action, scalar under thermodynamic 
coordinate transformations, which is stationary for general variations 
in the transport coefficients and the affine connection. 

A non-Riemannian geometry has been constructed out of the components of 
the affine connection, which has been determined by imposing the validity 
of the Universal Criterion of Evolution for non-equilibrium system relaxing 
towards a steady-state. Relaxation expresses an intrinsic physical property 
of a thermodynamic system. The affine connection, on the other hand, is 
an intrinsic property of geometry allowing to determine the equation for the 
shortest path. It is the author's opinion that a correct thermodynamical- 
geometrical theory should correlate these two properties. It is important to 
mention that the geometry of the thermodynamic space tends to be Rieman- 
nian for small values of the inverse of the entropy production. In this limit, 
we obtain again the same thermodynamic field equations found in ref. [3]. 
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The results established in refs [12|], for magnetically confined plasmas, and 



in refs 13|], for the nonlinear thermoelectric effect and the unimolecular 
triangular reaction, remain then valid. 

Finally, note that the transport equations may take even more general 
forms than Eq. (fT5j) . The fluxes and the forces can be defined locally as 
fields depending on space coordinates and time. The most general transport 
relation takes the form 

J M (r, t)= I dr' f dt'C^(X(r', t'))X v {v - r', t - t') (109) 
Jn Jo 

This type of nonlocal and non Markovian equation expresses the fact that 
the flux at a given point (r, t) could be influenced by the values of the forces 
in its spatial environment and by its history. Whenever the spatial and 
temporal ranges of influence are sufficiently small, the derealization and 
the retardation of the forces can be neglected under the integral 

C^ u (X{r',t'))X v {r-r',t-t') ~ £^(X(r,t))X»(r,t)5(r-r')5(t-t') (110) 

and the transport equations reduces to 

Mr,t)~T^(X(r,t))X v (r,t) (111) 



where 



(X(r, t)) = [ dr' ! dt'£^(X(r, t))S(r - r')S(t - t') (112) 
Jn Jo 



In the vast majority of cases studied at present in transport theory, it is as- 
sumed that the transport equations are of the form of Eq. (jllip . However, 
equations of the form Eq. (|109|) may be met when we deal with anomalous 
transport processes such as, for example, transport in turbulent tokamak 
plasmas [201] ■ Eq. (|110p establishes, in some sort, the limit of validity of 
the present approach: Eqs $63\) determine the nonlinear corrections to the 
linear ("Onsager") transport coefficients whenever the width of the nonlo- 
cal coefficients can be neglected. Moreover, it is implicitly understood that 
the thermodynamic quantities (temperature, pressure etc.) are evaluated 
making use of the local equilibrium principle. This last limitation, however, 
may be overcame by combining the TFT with the Extended Irreversible 



Thermodynamics (EIT) 21[. This will be addressed in the future. 
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6 Appendix 1: Transformation Law and Proper- 
ties of the Affine Connection Eq.(|55J). 

In this section we show that the affine connection Eq. (|55p transforms, under 
TCT, as in Eq. (i3Tj) and satisfies the postulates 1., 2. and 3. We first note 
that the quantity d^p+d^tpa transforms like a mixed thermodynamic tensor 
of third rank 

BX ,X BX P f)X u 



Thus, if Eq. (139j) transforms, under TCT, like Eq. (13ip . then so will be 
Eq. (|55p . Consider the symmetric processes. From Eq. (|22p . we have 



d 9'a8 dg pv dX? dXP dX u d 2 X p dX u d 2 X? dX u 



8X' K dXe dX' K dX' a dX'P + 9pu 8X' K dX' a dX'P + 9pu ~dX ,K dX'P dX 

(114) 

The thermodynamic Christoffel symbols transform then as 

Al' f r 1 dX' x dXP dX u dX' x d 2 XP 

avn avm av/R 



ap 1 pv J dX T dX /a dX'fi dXP dX' a dX'P 
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Recalling that a' = a, from Eq. (|114p we also find 



i 



-,X ,x O\g' c 



2a 



a/3) 



N Tr > dX' x dXP dX u 

—X n O{g pv ) dx/a —p 

'-x^J-^ 8 ^ 9 ^ 



(116) 



2a 



dX T dX ,a dX'P 



where Eqs (]20[) and Eqs (j27[) have been taken into account. Therefore, the 
affine connection 



pv 



T 

pv 

transforms as 



+ ^0{g pv 



r 



/A 

a/3 



ri, 



2(n+ l)a 
dX' x dXP dX v 



[8 T p X^O(g uv ) + 5lX r '0(g pu )] (117) 



dX' x d 2 XP 



pv dX T dX' a dX'P dXP dX' a dX'P 
Consider the general case. From Eq. (|22l) we obtain 



(118) 



l(d9L d 9'(3K d 9'a(3 

2\dX'P dX'P dX' K 



id 



dXe dXP dX u 
dX ,K dX' a dX'f 3 

d 2 XP dX v 

+9pv dX , adX ,/3 Q X ,K 



1VQ 



+ 



dg PQ dg, 



pv 



2\dXP dX v 



dXs 
(119) 



From Eq. (I24j) . we also have 

dfLu = df pv dx< dXp dxv 

dX'P dX< dX'P dX' a dX'P 
d fk = dfy, dX< dXP dX^ 
dX' a dXP dX'P dX' a dX'P 



d 2 XP dX r i 
+ ^QX'PdX' a dX'P 



d 2 XP dx** 
dX' a dX'P dX'P 



d 2 XP dX'v 
m dX'PdX'P dX' a 
d 2 XP dX r ' 



pv dX' a dX'P dX'l 3 
(120) 



Taking into account Eqs (|20j) and Eqs (|27j) we find 



X' K X'P 
X'X'p 



df 



dX'P 
dX' a 



xxri df P r,dXP dXP dX" 
Q dX v dX' K dX' a dX'l 3 

df uv dx? dXP dx» 

Q dXP dX'* dX' a dX'l 3 



+ X v XVf. 



+ x u xy, 



d 2 XP dX v 

PV Q X >a dX ,(3 Q X ,k 

d 2 XP dX v 



m dX' a dX'P dX' K 
(121) 



from which we obtain 



——X'X'vf 

2a' K \dX'P T 8X' a 



+ 



dXP dXP dX v 
dX' K dX' a dX'P 
1 _ d 2 XP 

+—X V X fprj — 



,L X X v[?Lpv 
2a 



■ d f»*)\ 
dX v dXP ) 



dX v 



a 



dX' a dX'P dX' 



(122) 
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where, as usual, a denotes the entropy production. Let us now re-consider 
the transformations of the following quantities 



dg' a » _ dg m dx< dxp dx^ 
dX'P dX< dX'P dX' a dX'P 
d 9pp _ dg, v dX< dXP 8X^ 
dX' a ~ dXP dX'P dX' a dX'v 



+9pv 



+9pv 



d 2 XP dX^ 
dX'PdX' a dX'P 

d 2 XP dX^ 
dX' a dX'P dX'v 



+9pv 



d 2 XP OX** 
dX'PdX'P dX' a 
d 2 XP dX^ 

dX' a dX'» dX'P 
(123) 



From these equations we obtain 



do' 

X^ +X 'P 



dX h 



xn Q_M_ , dg^\dx^_dxp_dx^_ d 2 xp 

dX< dXP ) dX'P dX' a dX'P p dX^dX' a 

(124) 

where Eqs (|27p have been taken into account. From Eq. (|124p we finally 
obtain 



1 
2a 1 



X K 



dg' , H 



+ 



ft 



dX'P dX' a 



1 

+-x p xy l 



PV 



dX h 



+ 



8XPJ 

d 2 XP 



dXe dXP dX u 



dX' K dX' a dX'P 
dX v 



a 



VI] 



dX' a dX'P dX 



Ik 



Summing Eq. ([119]) with Eqs (fl~22j) and ([125]) . it follows that 



VA 

a/3 



8X' X dXP dX v dX' x d 2 XP 



Vpv dX T dX' a dX'P 



+ 



dXP dX^dX'P 



where 



(125) 



(126) 



pv 



and 



N Te N, 



PQ 



5? 



pv 



N^x e X^df m df vv \ 
2a \dX v dXPJ 



2a 



dX v dXPJ 



with K 



PL> 



9 PS + —fprjX^Xg + —fp r] X v X p 



(128) 



Summing again Eq. (I127j) with Eq. (1113R and the first equation of Eq. (1116ft . 
we finally obtain 



r 



a/3 



T dX' x dXP dX u dX /X d 2 XP 
Tpu dX T dX' a dX'P + 
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dXP dX' a dX'P 



(129) 



where 

r T P » = r T pu + — x v o(g pu ) + s;^ + (130) 

It is not difficult to prove that the affine connection Eq. (|55p satisfies the 
postulates 1., 2. and 3. Indeed, if indicates a thermodynamic vector, we 
have 

A' x = A^^— (131) 



Deriving this equation, with respect to parameter we obtain 
dA' x _ dJP dX' x d 2 X' x dX T 

~dT ~ ~dT~dx^ + dx^dX'n & (132) 

Taking into account the following identities 

d 2 X' x dX' x dX' a d 2 XP dX' a dX'P dX' x d 2 XP 



dX T dX'n dXP dX T dX r idX' a dX T dX r > dXP dX' a dX'P 

(133) 

and Eq. (fT29l) . we find 

5A' X 5A^dX' x . 

~sT = ~JT~dx^ ( 134 ) 

The validity of postulates 2. and 3. is immediately verified, by direct com- 
putation, using Eqs (f3"2"j) and The validity of these postulates was 
shown above for a thermodynamic vector. By a closely analogous procedure 
it can be checked that the postulated 1. , 2. and 3. are satisfied for any 
thermodynamic tensor. 



Appendix 2: Derivation of the Thermodynamic 
Field Equations from the Action Principle. 



In this appendix, the thermodynamic field equations by the principle of the 
least action are derived. Let us rewrite Eq. ([62]) as 



/ 



R^sT - (r 



[IV 



f^)s 



Vgd n x 



(135) 



where the expression of S^ u is given by Eq. (|6ip . This action is stationary by 
varying independently the transport coefficients (i.e. by varying, separately, 
g^v and f^) and the affine connection T x . A variation with respect to T x 
reads 



Sir 



SRpug 1 



Vgd n x = o 



(136) 
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By direct computation, we can check that 

«V = OJT^V - (*r*,), A (137) 
Defining /C^ = ^/gg^ ', we have the identities 

(/c^^ A )|, = /cfr<5rJ A + /c^5rJ A|l/ 

(^<5r^)| A = ^5T^ + JC^JTA^ (138) 
Eq. (|136|) can be rewritten as 

J (K^ST^) lx d n X - J S^ST^d n X = (139) 
The thermodynamic covariant derivative of the metric tensor reads 

9aP\\ = 9a/3,A - TZ\9t,P ~ T }\3 V a ( 140 ) 

from which we find 

^ = -\g aP 9a P \x + \g ap g a p,x (ui) 

Taking into account that 5^/g = l/2^/ r gg fl ' / 5g IMU , Eq. (|14ip can also be 
brought into the form 

r^-^^ A + -Lv% = o (142) 

On the other hand, we can easily check the validity of the following identities 

(^^) |A = (x^<5rJj, A + (rj, - -L V5 )A + ^=y/grfivn% (143) 

Therefore, from Eq. (|142p . the terms 

| (/C^),^ and | (/C^<Sr*,)| A <fX (144) 

drop out when we integrate over all thermodynamic space. Eq. (|139p reduces 
then to 

SIv = ~J K^5T* x d n X + J JC^5T^d n X- J Sf 6T^d n X = (145) 
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It is seen that <5/r vanishes for general variation of 5T* U if, and only if, 

- \^x " + K$ - S?y/g = (146) 

Contracting indexes v with A, we find 

^:-Kp9 aP V~9 = ^ (147) 

where Eq. (I61j) has been taken into account. Thanks to Eq. (I147p . Eq. (|146() 

becomes 

= Kx9 Va V~9 + Kx9» a V9 (148) 
From the identity 5g^ u = —g fJ/a g ul3 5g a /3, we also have 

K\Z = y/9\ X 9T + Vgg?x u = \V99^ g aP g*p\x - Vgg^^g^x (149) 

Eq. (|148|) reads then 

- •W'W + \g aP 9ap\)^r = Kxg ua + Kxg^ (iso) 

Contracting this equation with we find, for 

rta/3|A = (151) 

where Eqs (|59|) have been taken into account. Eq. (|15U|) is simplified as 

- g^ a g up g a p\x = Kxg ua + Kxg^ (152) 

Contracting again Eq. (|152p with g m g vp , we finally obtain 

SV|A = -I^ASV - *°A5ar? (153) 

The first two field equations in Eqs (|63p are straightforwardly obtained con- 
sidering that from Eq. f)153|) we derive T^ v — T^ u = (see section [2|) . 
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